Testable polarization predictions for models of CMB isotropy anomalies 
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Anomalies in the large-scale CMB temperature sky measured by WMAP have been suggested as 
possible evidence for a violation of statistical isotropy on large scales. In any physical model for bro- 
ken isotropy, there are testable consequences for the CMB polarization field. We develop simulation 
tools for predicting the polarization field in models that break statistical isotropy locally through a 
modulation field. We study two different models; dipolar modulation, invoked to explain the asym- 
metry in power between northern and southern ecliptic hemispheres, and quadrupolar modulation, 
posited to explain the alignments between the quadrupole and octopole. For the dipolar case, we 
show that predictions for the correlation between the first 10 multipoles of the temperature and 
polarization fields can typically be tested at better than the 98% CL. For the quadrupolar case, 
we show that the polarization quadrupole and octopole should be moderately aligned. Such an 
alignment is a generic prediction of explanations which involve the temperature field at recombina- 
tion and thus discriminate against explanations involving foregrounds or local secondary anisotropy. 
Predicted correlations between temperature and polarization multipoles out to ^ = 5 provide tests 
at the ~ 99% CL or stronger for quadrupolar models that make the temperature alignment more 
than a few percent likely. As predictions of anomaly models, polarization statistics move beyond 
the a posteriori inferences that currently dominate the field. 



I. INTRODUCTION 

Several anomalies observed in the CMB data suggest 
a possible breaking of statistical isotropy on the largest 
angular scales. Broken isotropy would represent a radical 
revision of the standard cosmological model as statistical 
isotropy underlies all cosmological inferences. As such, it 
is important to test such claims by all means possible. 

Specifically, the measurements from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) [l|, [3, S 13 indicate 
an asymmetric distribution in power between the north- 
ern and southern ecliptic hemispheres Also, the data 
exhibit alignments in the temperature field; in particu- 
lar, the observed quadrupole and octopole of the CMB 
are planar and mutually aligned, and this plane is ori- 
ented roughly orthogonal to the CMB dipole direction 

m- 

The statistical significance of these preferred directions 
has been investigated by several authors [1, H, H, 
m, [3, Hi [Hj EBI J but the a posteriori choice of statis- 
tics makes their interpretation difficult [13, [3 ■ I'^" 
deed several generic tests of iso trop y y ield no evidence 
for statistical anisotropy [2^, [2l|, [U, [2^ . Since measure- 
ments of the temperature field have already reached the 
cosmic variance limit, further information from it will 
not be forthcoming, barring improvements in foreground 
and systematics modeling. On the other hand, large- 
scale CMB polarization measurements are still in their 
infancy and provide fertile grounds for future tests of 
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isotropy. While generic tests of isotropy have also been 
proposed for polarization [H, [2^ and unrelated statis- 
tical anomalies may be found, any physical model that 
purports to explain the temperature anomalies provides 
testable predictions for the statistics of the polarization 
field. Matching anomalies in polarization for such models 
can provide a means of going beyond both a posteriori 
inferences and blind searches for isotropy anomalies. 

In this Paper, we develop the radiative transfer tools 
required to investigate polarization statistics. Large- 
angle polarization arises from rescattering of CMB ra- 
diation from recombination during reionization. Given a 
model that predicts the statistical properties of the tem- 
perature field at recombination, Monte-Carlo realizations 
of the polarization field can be generated and analyzed 
for signatures of broken isotropy. 

We apply these techniques to models that break sta- 
tistical isotropy with a multiplicative modulation field 
[2^. These models posit that the observed temperature 
field is the product of a modulating field with superhori- 
zon wavelength fluctuations and a second field that car- 
ries fiuctuations at the observed scales. Gradients in the 
modulation field naturally pick out a preferred direction 
within our horizon. Such models have been proposed to 
explain both the hemisphere asymmetry 0, fl^ . [27| and 
the quadrupole-octopole alignment [26j . 

The structure of the paper is as follows. In § |lT] we 
discuss the general form of the modulation model and 
develop the simulation tools to predict the polarization 
field. The analytic framework for polarization predictions 
and explicit relations for the simplest modulation statis- 
tics are given in the Appendix. In ij lllll we study the po- 
larization predictions for the dipolar modulation models 
that are designed to explain the hemisphere asymmetry. 
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We continue in II VI by studying quadrupolar modulation 
models that explain the quadrupole-octopolc alignment. 
We discuss our results in § Ivl 



II. SPONTANEOUS ISOTROPY BREAKING 
A. Modulation Model 

We consider a model where a superhorizon scale modu- 
lation of the gravitational potential field causes the CMB 
temperature field to locally look statistically anisotropic, 
even though globally the model preserves statistical ho- 
mogeneity and isotrop y. We generalize the considera- 
tions of Gordon et al. j26| by explicitly working with a 
3D model of the gravitational potential as required for a 
polarization analysis. 

Let us begin by assuming that the source of gravita- 
tional potential or curvature fluctuations depends on the 
product of two fields (7(x) and /i(x). For example, these 
curvature fluctuations might arise from an inflaton de- 
cay rate that depended on two fields, each with their 
own quantum fluctuations. Let us further assume that 
the field /i(x) has only superhorizon scale fluctuations so 
that in any given Hubble volume, it takes on a determin- 
istic value whereas g{^) has subhorizon scale power that 
appears as stochastic fluctuations within the volume. 

Within our Hubble volume, an observer would sec bro- 
ken statistical homogeneity from the slow modulation of 
/i(x) across the volume. Furthermore, the local gradi- 
ent and curvature of /i(x) naturally picks out a direc- 
tion and breaks statistical isotropy as well. This bro- 
ken isotropy is then transferred to the CMB through 
the Sachs- Wolfe effect [2^] for the temperature field and 
Thomson scattering during reionization for the polariza- 
tion field. Nonetheless, a full spatial average over many 
Hubble volumes would show no statistical inhomogeneity 
or anisotropy. 

To make these considerations compatible with the 
observed statistical homogeneity and isotropy on small 
scales, we add another field that is responsible for fluc- 
tuations well below the horizon scale. Specifically, we 
model the Newtonian curvature fluctuation as 

$(x)=gi(x)[I + /^(x)] + <72(x), (I) 

where gi(x) and 52 (x) are Gaussian random fields and 
/i(x) is the modulating long-wavelength field. We will 
refer to the latter as the "modulating field" throughout 
this work. 

Several important properties of this model are revealed 
in Fourier space. The product of fields in physical space 
becomes a convolution in Fourier space 

<J>(k)=,9i(k) + j ^gi(k')Mk-k')+g2(k). (2) 

Two point correlations in gi and (72 arc determined by 



their respective power spectra 

(5r(k)5i(k')) - (27r)-^5(k-k')P,,(fc), 
(52*(k)52(k')) - (27r)35(k-k')P,,(fc). (3) 

We will typically assume that gi has power in fluctuations 
near the horizon scale and 52 has power on subhorizon 
scales. The modulating field h has only superhorizon 
modes. That the two point correlations do not couple 
modes of different k is a consequence of statistical ho- 
mogeneity and that the power spectrum is a function of 
fc = |k| only is a consequence of statistical isotropy. 

Notice that in the presence of the modulation ft,, $ does 
not obey statistical homogeneity and isotropy 

(<I>*(k)$(k')) = (2^)3<5(k-k')[P,,(fc) + P,,(fc)] 

+ [Pg,{k) + P,,{k')]h{\^ -k) (4) 

/li^k ~ 
^p,,(fc)r(k-k)Mk'-k), 

where the ensemble average is over realizations of gi and 
52 only. In particular, k modes that are separated by less 
than the horizon wavenumber are correlated and their 
amplitude depends on the direction of k despite the sta- 
tistical isotropy of gi. 

B. Temperature Field 

Broken homogeneity and isotropy are transferred onto 
the CMB temperature field through the Sachs- Wolfe ef- 
fect. Horizon scale fluctuations in the Newtonian cur- 
vature at recombination imprint large scale temperature 
perturbations as 

AT I 

— (x) = --<i>(x). (5) 

To the observer at the origin, temperature inhomo- 
gencities on a shell at the distance I?roc that a photon 
travels since recombination appear as anisotropics in a 
direction n 

AT AT 

— (n) = — (x- Accfi). (6) 

The analogue of a Fourier decomposition for the poten- 
tial field is the spherical harmonic decomposition of the 
temperature field 

AT ^ 

— (") = E E TemYUn), (7) 

e m=-e 

where Yem are spherical harmonics. 

Given a statistically homogeneous and isotropic poten- 
tial field, the two point correlations in the temperature 
field obey 

(^/'m'^*""') — Cf^Su'Smm' , (8) 



where Cj^ is the angular power spectrum of the temper- 
ature field. In the presence of the modulating field, the 
properties of the two point correlation of the potential 
in Eq. ^ carry over to the temperature field: modes of 
neighboring £ are correlated and the power depends on 
direction, i.e. the value of m. 

The two point correlation of the model in terms of the 
power spectra of the fields g and h was given in [2^ and 
is rederived in our framework of 3D fields in II In 
this framework, the modulation appears as the product 
of the Gaussian random field 51 (x) and the modulation 
field /i(x) to form the observed temperature field. For 
the temperature field, this is equivalent to the product of 
the 2D angular fields projected onto the recombination 
surface (see Fig. [1] solid circle) . Likewise the two point 
correlation can be expressed in terms of operations on the 
angular power spectra of the g and h fields in projection. 
We shall now see that these simplifications do not carry 
over to the polarization field. 



C. Polarization Field 

Broken homogeneity and isotropy are also transferred 
onto the CMB polarization field. A modulation model 
that is designed to explain temperature anomalies can be 
tested through the predictions it makes for polarization. 

Large-angle polarization is generated by the rescatter- 
ing of photons off free electrons during the rcionization 
epoch. A quadrupole anisotropy in the temperature of 
the radiation leads to linear polarization of the CMB. 
More specifically, from the temperature field at recombi- 
nation of Eq. ([5]) , the quadrupole moments at each posi- 
tion in space are given by 



AT 

T2™(x)= / dn'y;„(n') — (x'). 



(9) 



Note that T2m describes a set of 5 three dimensional 
source fields. Furthermore the angular integral is over 
the recombination surface of an (electron) observer at 
position X a distance D away from the true observer at 
the origin 




FIG. 1; Simulations of modulated temperature field, (a) 
Gaussian random unmodulated field g\. (b) Dipole modu- 
lated potential field with ^(x) = wi^/3/Atv{z/ Drcc) where the 
z axis is the modulation axis. Also shown (circle) is the sur- 
face at recombination upon which the observed temperature 
field originates. Note the power asymmetry in the modulated 
field. 



ADh' 



(10) 



where AD = Drcc — -D is the radius of recombination 
surface of the electron, as illustrated by Fig. [2l 

To the true observer at the origin, scattering at po- 
sition X = Dn generates a contribution to the Stokes 
parameters of the polarization in direction n given by 



-r{D) 



2 

X T2,n{Dh)±2Y2rn{n) , (11) 



m=-2 



where ±2Yim are the spin ±2 spherical harmonics and 
t{D) is the optical depth between a comoving distance 
D and the present epoch. When we specify r without 
the argument D, we refer to the total optical depth to 
rcionization. The analogue of temperature multipole mo- 
ments Tim for linear polarization are the Eim and Bim 
multipole moments 

(Q ± iU){n) = - ^ [Eem ± iS,,„] ±2l£,„(n) . (12) 

tni 

E and B modes are the tensor analogues of curl-frce and 
divergence free components of a vector. An E-tooAc has 
a polarization direction that is aligned with or orthog- 
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recombination 



FIG. 2: Schematic diagram for polarization generation. After 
reionization, free electrons along the line of sight n at a po- 
sition X see a quadrupole anisotropy on their recombination 
surfaces at a distance AD in their angular coordinates n'. 
Thomson scattering off these electrons produces the polariza- 
tion field viewed by the observer at the origin in direction 
n. 



onal to the direction that the mode amplitude changes. 
A B-mode has this direction rotated by ±45°. When 
polarization is generated by a scalar field, only i?-modes 
are formed. In Fourier space, this property is easily seen 
since the polarization direction from each k mode must 
follow the direction of k itself. Since the modulation 
preserves the scalar nature of the source of polarization, 
J5-modcs are absent even in the modulated case. Note 
that this is different from a direct modulation of the po- 
larization amplitude itself by a scalar quantity, e.g. from 
spatial fluctuations in r [29|. In that case, the polariza- 
tion direction is fixed by the potential field from recombi- 
nation but the amplitude carries additional directionality 
from the modulation. 

Predictions for the modulated polarization field, unlike 
the temperature field, require a 3D model of the modu- 
lation. Eq. ([To]) implies that the relevant modulation is 
at position x' in a sphere around each source at position 
X = Dn for each angular position n (sec Fig. [5]). Further- 
more, these sources lie in the interior of the recombina- 
tion surface D < Djoc of the temperature field. Likewise 
the polarization statistics cannot be expressed as an op- 
eration on the angular power spectra of the fundamental 
g and h fields corresponding to a single projection. The 
analytic framework for predicting the two point correla- 
tions is given in ^A 21 Given the cumbersome nature of 
this framework, we instead rely on numerical simulations 
to study modulated statistics as we shall now describe. 



We use the analytic framework to test these simulations 
(see ^A3p and to provide insight on statistical quantities 
that should be measured in the simulations and data (see 
El- 

D. Simulations 

Numerical simulations are conceptually the simplest 
and most general way of assessing the statistics of the 
modulated temperature and polarization fields. 

We first simulate the 3D potential field $(x) by com- 
bining realizations of the two Gaussian random g-fields 
with the deterministic modulation field h. We typically 
chose the length of the periodic box to be four times the 
diameter of the recombination surface with 512"^ pixels. 
This ratio of sizes and dynamic range allows us to cover 
the 2 < £ <10 scales of interest in this work. 

To obtain the temperature field, we then interpolate 
the field onto a 2D surface of radius Dree and weight it 
according to the Sachs- Wolfe relation Finally we use 
the HEALPix package [s^l to extract the multipole mo- 
ments of the field. The interpolation is done in every unit 
vector n corresponding to 12, 288 HEALPIX pixel centers 
at resolution A'sidc = 32. As an example, slices through 
simulations of unmodulated and dipole-modulated tem- 
perature fields are shown in Fig.[TJ The statistics of these 
temperature simulations match those of the 2D approach 
where the modulation is directly in angular space on the 
recombination surface. We will occasionally use 2D sim- 
ulations where only the temperature field is required and 
for unmodulated cases to build large samples. 

The polarization field is equally straightforward to sim- 
ulate but numerically more involved. We take the same 
3D simulations employed for the temperature field and 
compute the quadrupolar temperature sources for 80 red- 
shifts in equal steps along each line of sight, and in every 
unit vector n corresponding to 12, 288 HEALPIX pixel 
centers at resolution A'sidc ~ 32 around the sky. Then, we 
generate maps of polarization by means of Eq. (jlip . Once 
the maps are created, the coefficients Eim arc obtained 
using the HEALPix package [l^l, and various polariza- 
tion statistics described below are computed. 

We test this pipeline against analytic calculations of 
unmodulated and simple modulation models in f|X3l It 
is the main tool used to assess the temperature and po- 
larization statistics in this work. It can be used for arbi- 
trarily complicated modulation fields and any statistical 
property of the temperature and polarization fields. 



III. DIPOLAR MODULATION 

Eriksen et al. Q showed that there is a hemispheri- 
cal or dipole asymmetry in the amount of power in the 
WMAP data at muhipoles £ - 20 - 40. The axis of this 
asymmetry points in a direction that is close to the South 
Ecliptic Pole (SEP) making the power in the southern 
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hemisphere larger than the north. It was subsequently 
pointed out that a dipolar modulation of the form of 
Eq. (|Aip with a modulation field of 



h{x) = wiYj 



10 



(13) 



could explain the asymmetry though not the curi- 
ous ahgnment with the SEP 0, [H, Note that 

Yio(^/-Drcc) = ^3/47r(z/I?roc) in a coordinate system 
where the z-axis points along the modulation axis. We 
set the g2 field to have power only at £ > 40 and hence its 
contribution is irrelevant in our considerations. We take 
the gi field to have a scale invariant power spectrum 



k'PM 9 



27r2 



25 



(14) 



with (5^ = 4.6 X 10~^ which corresponds to the WMAP 3- 
year normalization for a cosmology with matter density 
ri„i = 0.24, cosmological constant 51a = 0.76, baryon 
density flbh"^ = 0.022, Hubble constant Hq = 73 km 
Mpc~^, and spatially flat spatial geometry. 

Once wi is specified, we can use our simulation pipeline 
to make realizations of the polarization field. An extreme 
example with wi = 2.5 is shown in Fig. [3] Note that 
the hemispherical power asymmetry carries through to 
the polarization. The WMAP temperature field however 



in not only TT 

{Tg„,Te+i,,n) oc wi , (15) 
but also in all of the two point combinations involving E 



(16) 



specifically XY e {TE, ET, EE}. We derive the propor- 
tionality coefficients in Appendix lA 41 Each multipole 
pair then forms a noisy estimator of the dipole modula- 
tion parameter wi. Even in the absence of detector noise 
and foregrounds, the cosmic variance of the temperature 
and polarization fields provide noise to the measurement. 
Eqs. (|A39p and (|A4G[) give the unbiased minimum vari- 
ance combinations of the multipole pairs as estimators of 

Wi. 

The first question we ask is: if the observed dipole 
asymmetry at € ~ 20 — 40 is a statistical fluke and that 
the true ensemble average is wi = 0, how well can tem- 
perature and polarization at € < 10 rule out a finite value 
of wi near the best fit of wi « 0.2? Under this assump- 
tion, we show the distribution of the XY estimators and 
the joint estimator in Fig. |4l These distributions were 
generated from 10^ 2D simulations of the unmodulated 
temperature and polarization fields in accordance with 
the power spectra of Eq. (K9\i . (|A20|) and (|A30|) . 




FIG. 3: Map composed of the first 10 multipoles of the E 
polarization field in a dipole modulation model with wi = 2.5. 
The power is larger in the southern hemisphere. 
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suggests a much weaker dipole modulation that can only 
be detected statistically in the polarization. We will base 
our parameter choices on the analysis of Ref. |12] which 
found wi = 0.23 ± 0.07 with wi ^ excluded at the 99% 
CL. For illustrative purposes we will take a reionization 
model with t = 0.09 where hydrogen is fully ionized out 
to a maximum redshift and neutral thereafter. We test 
sensitivity to the reionization history in ^IIIBI 



FIG. 4: Probability distribution of the dipole modulation es- 
timators for a model with uii = where XY represents 
the two fields whose £ and £+1 multipoles are correlated: TT 
(red long-dashed curve), TE (black dotted-dashed curve), ET 
(violet dotted curve) and EE (blue thick dashed curve). The 
minimum variance joint estimator wi is shown as the solid 
curve. We have assumed cosmic variance limited measure- 
ments and ^max = 11 for all estimators. 



A. Dipole Modulation Statistics 

The dipolar modulation in the potential field gives rise 
to couplings between £ and £' = £ ± 1 multipole moments 



Given TT alone, the hypothesis that wi ^ 0.2 can 
only be weakly tested. Specifically, a value of wi > 0.2 
can be tested at the 86.4% confidence level (CL) by TT. 
Each of the polarization estimators can in principle test 
dipole modulation at a comparable level: wi > 0.2 can be 
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tested at 81.6% confidence by EE, at 83.4% CL by TE 
and at 82.1% CL by ET. The combination of the four 
estimators however is more powerful: wi > 0.2 can be 
tested at 98.4% CL by the joint estimator. For reference, 
the central value of [I4I, wi = 0.23, yields a 99.3% CL 
test. Furthermore the 95% upper limit on wi achievable 
with the joint estimator is wi < 0.15. 
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FIG. 5: Standard deviation of the individual estimators Wi 
and the joint estimator (black) as a function of maximum 
multipole pair employed in the estimators ^max- The indi- 
vidual XY pairs are TT (red long-dashed curve), TE (black 
dotted-dashed curve) , ET (violet dotted curve) and EE (blue 
dashed curve) from the analytic calculation of the Appendix. 



Most of the polarization information comes from mul- 
tipoles in the range of ~ 4—8. In Fig. Owe show the rms 
fluctuation in the estimators as a function of ^max, the 
maximum multipole employed in the estimators. Note 
that the estimators involving polarization saturate by 
10. 

The next question we ask is: given a true model of 
wi = 0.2, how well can an isotropic model of wi = 
be excluded? Fig. [5] shows the distribution of wi com- 
ing from 545 realizations of a 3D modulated model (his- 
togram, ^max = 11)- In these simulations wi < oc- 
curred in less than 1% of the realizations. Furthermore 
the mean of the samples (wi) = 0.20 is consistent with 
being unbiased and the standard deviation a = 0.08 is 
consistent with the expectations from Fig. [5l 

B. Model Dependence and Noise 

These tests are only weakly sensitive to the assumed 
reionization history for a cosmic variance limited observa- 
tion. The main difference is that a linear modulation pro- 
duces a weaker effect for nearby sources of polarization 
and hence the signal decreases with r. As an example, 
we calculate the level of detectability for an unphysically 
small optical depth of r = 0.0279. Even for this extreme 
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FIG. 6: Histogram of the joint estimator wi coming from 545 
3D realization of the dipole modulation model with wi = 0.2 
and T = 0.09. For comparison, the probability distribution 
for the isotropic wi = model from Fig. |4] is also plotted 
(smooth curve). 



case, values of wi > 0.2 are ruled out at the 96.4% CL by 
the joint estimator. An intermediate case of t = 0.055 
rules out values of wi > 0.2 at the 98% CL. 

Because most of the signal comes from low multipoles 
in the polarization £ 4 — 8, the cosmic variance limit 
described in the previous section is within the statistical 
reach of the Planck satellite. For definiteness we take 
the 143 GHz and 217 GHz frequency channels and com- 
bine them into an effective white noise power spectrum 
of Nf^ = (67.6/iK-arcmin)^ at low multipoles. Noise in 
the temperature field is negligible compared with cosmic 
variance. Quantitatively, for a reionization history with 
optical depth t = 0.09, and the noise level considered 
for Planck, values of wi > 0.2 are ruled out at 97.8% 
CL by the joint estimator wi. A space-based experi- 
ment with no foregrounds, 5 frequency channels and a 
noise level of 18/^K-arcmin per frequency channel would 
essentially saturate the cosmic variance bound in each 
channel. Given these results, foregrounds and systemat- 
ics will probably be the limiting factor in how well dipole 
modulation can be tested by polarization. We leave this 
subject to a future work. 

Finally, the model spectrum assumes that the tempera- 
ture anisotropy comes purely from the Sachs- Wolfe effect. 
A ACDM model would predict that the ensemble aver- 
aged quadrupole has a substantial contribution from the 
so-called Integrated Sachs- Wolfe (ISW) effect from the 
decay of potentials during the acceleration epoch. Under 
ACDM, the measured low power quadrupole anisotropy 
could indicate that either the ISW contributions to the 
quadrupole are by chance small on our sky or that there is 
a chance cancellation between the two contributions. Fi- 
nally, it could indicate an alternate model of the dark en- 
ergy where the ISW contributions arc smaller [3l|, [H, [s^ 
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or designed to cancel the Sachs- Wolfe effect [34|. In all 
of these explanations, the polarization field is unaffected. 

Nonetheless, the dipole modulation tests involving the 
temperature field are slightly weakened if the ISW con- 
tributions to ACDM are added back as a source of noise. 
Correspondingly the joint estimators also weaken. For 
example, for a reionization model with optical depth 
r = 0.09 the joint estimator can test values of wi > 0.2 
at the 97.6% CL. 

In summary, prospects for future polarization tests of 
dipole modulation at the wi w 0.2 level at the ^98% CL 
are fairly robust to model assumptions. 



IV. QUADRUPOLAR MODULATION 

de Oliveira-Costa et al. Q pointed out that there 
arc several striking features of the temperature dipole, 
quadrupole and octopolc in the WMAP data. The oc- 
topole is highly planar and defines a preferred axis that 
is near the dipole axis. Moreover the quadrupole is also 
weakly planar in the same coordinate system. 

Gordon et al. [2^ subsequently showed that a 
quadrupolar modulation with a preferred axis that is by 
chance near the dipole direction can make these features 
substantially more likely. Specifically, we will take a mod- 
ulation to the potential field of the form of Eq. ^ , with 



(17) 



Note that the z-axis required here is not the same as 
in the dipole modulation case. Since Y2q{z / D-^cc) = 
(l/4)-\/5/7r [3(z/Z?rcc)^ — l] , this introduces a quadratic 
term to the spatial modulation. 

These anomalies are confined to the low multipoles and 
so the modulation must also be confined to low multi- 
poles. For the angular modulation considered in [26j . 
this was simply imposed as a sharp suppression in mul- 
tipole space of the gi field. Spatial modulation is more 
complicated in that even a sharp suppression in Fourier 
space is broadened by projection effects. Given these 
considerations, we take a two field model defined by 



27r2 

k'^ Pg^ (fc) 

27r2 
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ki<k< k2 , (18) 

k > fcg , (19) 



where the normalization S,^ is given as in the dipole mod- 
ulation model by the WMAP normalization. We intro- 
duce a scaling parameter / that controls the amount of 
power at low k relative to this normalization. We begin 
our analysis by finding regions in the parameter space 
of {^2, fci, fc2, fcs, /} that maximize the probability of the 
observed quadrupole-octopole alignment. 



A. Maximizing Temperature Alignments 

To quantify the quadrupole-octopole alignment, de 
Oliveira-Costa et al. Q introduced the normalized an- 
gular momentum along a given axis {9, cj)) 



(20) 



where the T^'^ are the multipolc moments of the tem- 
perature field in that preferred frame. More specifically, 
from the moments defined in galactic coordinates Tim, 
the rotated ones are 



T' 

-'-hn' 



Tim 



(21) 



where the Euler angles are ip G [0, 27r], 9 G [0,7r], and 
(f) S [0, 27r] and D^^,{~'iIj, ~9, —(f)) is the Wigner rotation 
matrix [2^ . Note that Lf is invariant under the final 
rotation by ip. 

In particular, the statistic 
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captures both the alignment between the quadrupole and 
octopole, and the planarity of the octopole. The pre- 
ferred axis is then chosen as the one which maximizes this 
statistic. For the ILC map of the 3-year WMAP data cor- 
rected for the kinematic quadrupole, = 0.943. This 
is slightly smaller than for the 1-year WMAP data in 
the TOH map used in previous studies: L23 = 0.96 
BUEl- For an unmodulated model, £23 > 0.943 occurred 
in only 0.85% of 10** 2D realizations. 

Quadrupole modulation can dramatically increase the 
probability of a high angular momentum statistic. To 
find the parameters that maximize this probability, we 
first vary the parameters {w2, ^1,^:2} with fca — > 00 and 
find the set that gives maximum number of realizations 
with 1/23 ^ 0.943. Taking k^ 00 eliminates the con- 
tributions from the unmodulated (72 field that could oth- 
erwise destroy alignments at low multipoles. The align- 
ment probability then is also independent of /. 

We then generate 10^ 2D Monte Carlo (MC) reahza- 
tions for each of the modulated models. For each real- 
ization, the estimator L23 is computed on a grid of 10® 
positions on the celestial sphere by performing the rota- 
tion of Ti„i in Eq. (|21[1 . Taking the maximum value of 
1/23 attained for each realization, we compute for each 
model the percentage of realizations with L23 > 0.943. 

The region of parameter space that corresponds to the 
highest alignments is centered around W2 ~ —7, in agree- 
ment with [2Q] . Furthermore, wavenumbers that best 
project onto a multipole i correspond to fcl?rec ~ ^; more 
specifically fcDroc ~ 4.2 for the region that best projects 
onto the quadrupole and octopole at the distance to re- 
combination. Furthermore keeping ki and k2 in a narrow 
range around this value eliminates contamination from 
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the modulation of longer and shorter wavelengths in the 
quadrupole and octopole. 

For definiteness, we take W2 — —7, kiD^cc ~ 4.2 and 
fe-Dioc = 4.23 though any {k2~ki)D^cc ^ 0.05 would pro- 
duce similar results due to unavoidable projection effects. 
With this choice, 21.8% of reahzations have > 0.943. 
This number is smaller than that quoted in |26| due to 
these projection effects but is still a factor of 25.6 larger 
than the probability in unmodulated models. 

While taking ks —^ oo maximizes alignments, it does 
not provide for a realistic model. The absence of high 
wavenumbers would be in conflict with the observed fluc- 
tuations at higher multipole moments. This is exhibited 
in a poor ILC likelihood 



/:(T|c) = 



1 



(27r)^/2Vd^ 



exp[--TtC-iT] 



(23) 



where T is the N element vector of observed T^,„ values 
and C is their signal plus noise covariance matrix. To 
remedy this problem we keep the first field fixed to the 
parameters above but add in the second field at a finite k^ 
and maximize the likelihood over their relative amplitude 
/■ 
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FIG. 7: Effect of the unmodulated field on alignments and 
goodness of fit of quadrupole modulation models. The upper 
panel shows the fraction of realizations that have a higher 
angular momentum statistic L23 than WMAP 3-year ILC as 
a function of fcaDrec which controls the contribution of the 
unmodulated field to the low multipoles. The dependence of 
Axeff (relative to ACDM, dashed line) is shown in the lower 
panel. 



Fig. [7] (lower panel) shows the dependence of XcS = 
-2 In/: on fcaDi-cc for £ = 2 - 5 relative to ACDM (see 
below Eq. ^^). Decreasing this parameter creates a 
better fit to the data at the expense of decreasing the 
probability of the quadrupole-octopole alignment (upper 
panel) . As a compromise we choose a model with Ax^g = 
that neither improves nor worsens the fit. The model 



parameters then become 

fciAec=4.2, fcai'rcc = 4.23 , 

W2 = ~7, / = 5.77, r: 



fcsAec = 6.697, 
0.09, (24) 



for which the probability of L23 > 0.943 is 16.6% (see 
Fig.©. 
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FIG. 8: Distribution of L23 (red solid line) for the quadrupolar 
modulation model of Eq. (|24|) . 16.6% of the 10* realizations 
of this model have L23 > 0.943. The distribution for the 
isotropic case is shown as a dashed line for comparison, where 
0.85% have Lig larger than the WMAP ILC value: 0.943 
(vertical line). 



Fig. [HI shows the temperature angular power spectrum 
of this model along with the WMAP 3-year temperature 
data for the first 10 multipoles. The bands shown here 
represent the 68% and 95% cosmic variance confidence 
regions for this model taken from 10"*^ 2D realizations of 
the temperature field. Note that the cosmic variance is 
larger for the modulated model due to the m-dependence 
of the power. 

In Fig. [TO] we show a realization of the model that 
has a high normalized angular momentum L23 = 0.94, 
shown in the frame of the modulation, illustrating how 
the modulation takes a random isotropic sky and aligns 
the quadrupole and the octopole. 



B. Polarization Angular Momentum 

With the quadrupolar modulation model set by 
Eq. (I24p to maximize temperature alignments and fit the 
WMAP temperature data, we can now use our simula- 
tions to make polarization predictions. 

We first examine the normalized angular momentum 
statistic of the polarization field. We define this statistic 
in the same way as for the temperature field 



f 2(pol) 



(25) 
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FIG. 9: Temperature angular power spectrum of the 
quadrupolar modulation model of Eq. (|24|l (solid line). The 
WMAP ILC data with noise errors are also shown. The solid 
bands correspond to the 68% and 95% cosmic variance con- 
fidence regions of the model. The dashed line represents the 
ACDM model against which this model is compared in Fig. [7] 



OUADRUPOLE OCTOPOLE 

FIG. 10: Quadrupole and octopole of the temperature field 
from one realization of the quadrupole modulation model of 
Eq. (|24|l . The normalized angular momentum here is Lis = 
0.94. The range goes from — 50/^K to 50/iK for both cases. 



with the combined quadrupole and octopole statistic be- 
ing 

^1^°'^ = ^(^2^''°'^+L^(P°')). (26) 

There is one important difference between this statistic 
and the temperature based one. In the case of polariza- 
tion predictions, the temperature field has already de- 
fined the preferred axis against which we should measure 
the polarization angular momentum. Consequently, no 
search over directions is required. 

Fig. [TT] shows a histogram of L2^P°'^ for 571 out of 
3482 3D realizations of the modulated model that have 
as high or higher temperature alignment than WMAP 
L^g > 0.943. We first compare this with the same statis- 
tic for 10^ 2D isotropic realizations of an unmodulated 
polarization field (red dotted line). There is a tendency 
for the polarization field to be aligned compared with 
random isotropic realizations. An example of a typical 
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FIG. 11: Histogram of -£23^°''' for 571 3D realizations of 
the quadrupole modulation model of Eq. (|24|l that also have 
-£23 ^ 0.943. For comparison we also show the distribu- 
tion for isotropic random polarization fields (red dotted lines). 
Isotropic realizations constrained to the Tim values of the ILC 
map have a distribution shifted to higher alignment values 
given the temperature-polarization correlation with just the 
Sachs- Wolfe effect (blue dashed lines) and with ACDM ISW 
contributions added (black solid curve). 




OUADRUPOLE OCTOPOLE 

FIG. 12: Realizations of the quadrupole (left column) and 
octopole (right column) of the E component of the polariza- 
tion field in order of increasing alignment: isotropic case with 
l2(p°') = 0.46 (first row), quadrupolar modulation model with 
parameters given by Eq. (|24|l with L^^"^^ = 0.67 (second row) 
and an isotropic Sachs- Wolfe model constrained to match the 
WMAP ILC temperature muhipoles with -£23""°'' = 0.9 (third 
row). The color scale goes from — 0.3/iK to 0.3/^K in all pan- 
els. 
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realization is shown in Fig. [12] (second row) and is to be 
compared with a typical isotropic realization (first row). 
Note the similarity between temperature alignments and 
E polarization aligments in Figs. [TUl and [T^ Both have a 
deficit of power near the poles in the frame of the modula- 
tion. However the polarization amplitude {Q^ + U^y^^ is 
in fact the largest at the poles. In this direction, the elec- 
tron scattering sees the planar quadrupolc in the trans- 
verse plane and scatters it into a maximal polarization. 
Foregrounds which have high polarization where there is 
high intensity will not have this pattern. 

However, the expected correlation between the temper- 
ature and polarization fields implies that there should be 
some tendency for polarization alignment in any expla- 
nation in which the alignment of the temperature field re- 
flects contributions from the recombination surface, and 
not a foreground, secondary or systematic effect. This 
includes an isotropic unmodulated model in which the 
temperature alignment is a chance occurrence. 

To quantify this tendency, we consider two isotropic 
models which have been conditioned to have = 0.943 
for the temperature field. More specifically, we draw 
a polarization field that is consistent with the observed 
temperature field Tim and model power spectra. 



a 



TE 



(jTT 



Tp, 




EE 



r<TT 



(27) 



where n^m is a complex Gaussian variate of zero mean 
with (nl^ntm) = 1 and n}^ = (-l)'"n^,_m. 

The first is the Sachs- Wolfe model where all of the tem- 
perature anisotropy arises from recombination. Under 
this assumption, the temperature alignments are a real 
feature of the recombination surface and the ISW effect in 
our local universe happens to be small in the quadrupole 
and octopole. Fig. [11] (blue dashed lines) shows the po- 
larization statistic for 10^ 2D constrained realizations of 
this model. An example of a highly aligned realization is 
shown in Fig. [T^] (third row). 

Compared with the modulated model, there is an even 
higher probability of polarization alignment here. The 
reason is that in the former, the strength of the spa- 
tial modulation increases as the square of the distance. 
The closer distance to rcionization implies a substan- 
tially weaker modulation of the polarization field as com- 
pared with the temperature field. Ironically, a very 
high value of L^'"^"^^ would strongly disfavor the mod- 
ulated model compared with the isotropic Sachs- Wolfe 
model. The 95% upper limit for the modulated model is 



f 2(pol) 
^23 



0.77. 



It is also possible that the alignment of the quadrupole 
and octopole does not represent an aligned spatial con- 
figuration at all and is instead the result of a chance 
combination of the Sachs- Wolfe effect and ISW effects. 
These effects arise from very different physical scales but 
are seen in projection at the same angular scale. To il- 
lustrate this case, we draw 10^ 2D realizations of the 
ACDM model (see below Eq. that are similarly con- 



ditioned to have I/23 = 0.943 for the temperature field. 
The histogram of the polarization statistic is shown in 
Fig. [11] (black sohd curve). The presence of a tempera- 
ture contribution that does not correlate with polariza- 
tion shifts the distribution downwards and decreases the 
probability of polarization alignments. The distribution 
is still distinct from the truly random isotropic simula- 
tions (red dotted lines). The 95% lower limit on this 



model is E 



2(pol) 
23 



> 0.39. 



The final possibility is that the temperature 
quadrupole and temperature field is nearly all local sec- 
ondary anisotropy and that the Sachs- Wolfe effect is 
anomalously low (e.g. [3^, H^). The polarization dis- 
tribution would then look like the truly random isotropic 
simulations (see Fig.[TTl red dotted line). In this case, the 
polarization amplitude itself would also be anomalously 
low. 

In summary, a very low value of E^'^'"^^ ^ 0.4 with 
a normal polarization amplitude would suggest that the 
origin of the temperature alignments is not cosmological. 
A very low value with suppressed polarization amplitude 
would imply a purely secondary origin to the temperature 
alignment. A very high value L^''^"^^ ^0.8 would rule out 
the modulated model and also imply that the ISW effect 
is nearly absent in the WMAP low multipoles. 
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FIG. 13: EE and TE power spectra of the modulated model 
of Eq. p4]l (red dashed lines) along with the ACDM power 
spectrum (black solid lines). The modulated model has a 
lower polarization at the same r(= 0.09). The units are 



C. Quadrupole Modulation Statistics 

Beyond the angular momentum of the quadrupole and 
octopole, a quadrupolar modulation model makes further 
testable predictions. These predictions however are more 
sensitive to the choice of model parameters. 
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First, the angular power spectrum of the polarization 
and temperature cross correlation differs in its relation- 
ship to the temperature power spectrum due to the mod- 
ulation. In Fig. [ini we show these spectra for the model 
of Eq. (|24p . Note that for the same optical depth r, the 
polarization spectra have less power. This is a conse- 
quence of the decrease in modulation amplitude at the 
closer distance to reionization. Unfortunately, this sig- 
nature is largely degenerate with a change in the optical 
depth and so is not a unique signature of modulation. 

Similar to dipole modulation, quadrupole modulation 
generates a more unique signature in the correlation be- 
tween neighboring multipole moments 



(28) 



where XY £ {TT, TE, ET, EE). In our two field model, 
most of the modulation is confined to the low multipolcs 
and so the signal is mainly in + 2 = 4, 5. Note that an 
isotropic model that by chance has quadrupole-octopole 
temperature alignments will still predict a zero expecta- 
tion value for these quantities. 

Let us define a correlation statistic 



.771 I 



(29) 



which tests the phase alignment of the moments and 
takes on values from — 1 < sf^ < 1. In the absence 
of modulation, these statistics have zero mean but any 
given pair has large cosmic variance. Given a modulation 
model that predicts 



IZXY\ 
\°e.m /mod ; 



we can construct a weighted sum 



,XY 



\ ^m,mod j 



(30) 



(31) 



where ^ = 2, 3. This statistic has the property that its 
expectation value is zero for an isotropic field and 1 for 
the chosen modulation model. 

We begin by considering the XY = TT case and the 
modulated model of Eq. ((24|) as in the previous section. 
The distribution of s is shown in Fig. [T3] (dashed lines, 
right curve) from 10^ 2D simulations. The WMAP ILC 
data has s = 0.076 and disfavors the modulated model 
but not at high significance: 3.5% of the modulated sim- 
ulations have s less than the data. Conversely the data 
would be a typical realization of an isotropic model (see 
Fig. [Ml dashed lines, left curve). Likewise an s < 
occurred in 2.5% of the modulated simulations, and an 
s > 1 occurred in 2.3% of the isotropic simulations. 
These two tails of the respective distributions represent 
the power of the statistic to distinguish typical realiza- 
tions of the models. 

The addition of polarization information would make 
these tests definitive by reducing the probability in the 



tails (see Fig. [Mj solid curves). When all fields are con- 
sidered XY e {TT,EE,TE,ET}, the probabilities in 
the tails drop to s < in only 0.09% of the 3482 modu- 
lated 3D simulations and s > 1 in 0.24% of the 10^ 2D 
isotropic simulations. The latter probability assumes a 
Sachs- Wolfe only model for the temperature field. If the 
ISW contribution of ACDM is factored in, the probability 
increases slightly to 0.34%. 

Consequently a quadrupolar modulation with parame- 
ters given by Eq. ([M]) will be definitively tested through 
the s statistic by polarization measurements that ap- 
proach the cosmic variance limit for i < 5. However a 
model with a lesser amount of modulation at £ = 4, 5 
may evade such constraints. The amount of modula- 
tion in these multipolcs is controlled by and hence s 
mainly constrains this parameter. In the model we have 
been considering, k^Drcc = 6.697. Lowering this value 
increases the contribution of the unmodulated field for 
£ < 5. But the price is to make the quadrupole-octopole 
alignment less likely (see Fig. [71). 

As a specific example, for k^Dj-oc = 5.16 and / = 5.13, 
the probability of L23 > 0.943 decreases from 16.6% to 
10% (see Fig.[7l). On the other hand, it improves the fit 
to the WMAP data by making Axlg = — 9 relative to 
ACDM. Correspondingly, the s-statistic also weakens in 
its discriminating power. For the TT case, the WMAP 
ILC data has s = 0.065. For this case, 10.6% of 2521 
simulations have a smaller value of s. Similarly, a value 
of s < occurred in 9.1% of the modulated simulations, 
and s > 1 occurred in 10.6% of the isotropic simulations. 
With the addition of polarization information these tests 
become more significant: s < only occurred in only 
0.24% of the modulated simulations and s > 1 only oc- 
curred in in 0.50% of the isotropic simulations. 

In summary, while the s-statistic cannot rule out 
quadrupolar modulation in general, it can severely limit 
the amount of modulation allowable at ^ < 5 and hence 
test its ability to explain the quadrupole-octopole align- 
ment. 



V. DISCUSSION 

Statistical anomalies in the large-scale CMB temper- 
ature sky measured by WMAP have been suggested as 
possible evidence for the violation of statistical isotropy 
on large scales. In any physical model for such a viola- 
tion, there are testable consequences for the CMB polar- 
ization field. 

We have developed a radiative-transfer simulation- 
based approach to predicting these consequences. We 
have applied this approach to modulation models of bro- 
ken isotropy. In this class of models, a long-wavelength 
modulation of shorter-wavelength fluctuations breaks 
isotropy locally, while the universe remains globally 
isotropic and homogeneous. 

A linear or dipolar modulation can be used to explain 
the observed power asymmetry between the northern and 
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FIG. 14: Distribution of ^±2 correlation statistic s. Distribu- 
tions on the right with {§) — 1 represent the modulated model 
of Eq. (|24p for the joint estimator (black solid lines) and the 
temperature estimator (red dashed lines). Curves on the left 
with (s) = represent the isotropic Sachs- Wolfe model in the 
same manner. The histogram comes from 3482 3D realiza- 
tions of the modulated model and other distributions from 
10^ 2D realizations. The vertical red dashed line represents 
the observed value for temperature only for this modulated 
model, using the WMAP ILC map: s = 0.076. 



southern ecliptic hemispheres at multipolc moments £ ~ 
20 — 40 [3, 0, m, [131 . Linear modulation produces a 
dipole coupling that correlates multipole moments £ to 
^ ± 1 in all of the two point functions of the temperature 
and E polarization field. Using an analytic framework 
developed in the Appendix, we derive an estimator of 
the modulation amplitude from these fields from these 
correlations. 

We find that the addition of polarization information 
substantially improves the ability to test for dipole mod- 
ulation at the lowest multipoles £ < 10. Under a wide 
range of assumptions, the dipole modulation explanation 
of the hemisphere asymmetry can be tested at the 98% 
CL or greater. These tests will require near cosmic vari- 
ance limited polarization measurements at the low multi- 
poles. Such measurements are within the statistical reach 
of the Planck satellite but systematic contamination from 
foregrounds will have to be sufficiently controlled. 

A quadratic or quadrupolar modulation in the temper- 
ature field can be used to explain the alignment of the 
quadrupole and octopole [1, [2^ . We generalize the con- 
siderations of Ref. [20] and develop full 3D models that 
enhance the probability of alignments by a factor of ~ 20. 

These models predict that the quadrupole and oc- 
topole of the i?-polarization will tend to align as well. 
However alignment of the polarization mainly tests that 
the temperature alignment is cosmological and originates 
from recombination, not the specifics of the modulation 
model. As such its main use will be to reject explana- 



tions of the alignments due to foregrounds or purely lo- 
cal secondary anisotropy. In fact a very high alignment, 
comparable to the temperature field, is not expected in 
the modulation model due to the smaller strength of the 
modulation at the distance to reionization. If such an 
alignment is seen then it is more likely that it arises from 
a chance configuration of the temperature field at recom- 
bination. 

Quadrupolar modulation models can be distinguished 
in a manner similar to the dipole modulation but through 
a correlation between multipoles £ and ± 2 in the tem- 
perature and polarization fields. We find that the combi- 
nation of temperature and polarization information out 
to £ < 5 can strongly constrain the level of modulation in 
the low multipoles. As a consequence, polarization mea- 
surements can in principle rule out models that make the 
temperature alignments more than a few percent likely 
~ 99% CL or stronger. 

In general, matching signatures of broken isotropy in 
the polarization provide a powerful test of competing ex- 
planations of the temperature anomalies. Explanations 
that involve foregrounds (e.g. [s^]) or local gravitational 
effects H^lll] would tend to predict different signatures. 
This is especially true since matching anomalies in the 
E polarization are non-local in the Stokes parameters. 
Possible experimental systematics would be specific to a 
given experiment, and would thus have to be analyzed in 
a case by case basis. Therefore we do not consider them 
in this paper. 

The examples we have studied illustrate that polariza- 
tion in the future will provide a fruitful testing ground 
for explanations of the observed temperature anomalies. 
As predictions of anomaly models, polarization statistics 
move beyond the a posteriori inferences that currently 
dominate the field. 
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APPENDIX A: ANALYTIC MODULATION 
FRAMEWORK 

In this Appendix, we construct an analytic frame- 
work for modulated temperature ( ^A 1[) and polarization 
statistics (f|X2|) and evaluate them explicitly for dipole 
modulation. We use these analytic calculations to cross 
check numerical results in 31 and to construct estima- 
tors of the dipole modulation in ^A4\ 



1. Modulated Temperature Field 

For the temperature field, the effect of modulation can 
be viewed as either the direct projection of the modulated 
(convolved in Fourier space) field or the convolution in 
harmonic space of the multipole moments of the origi- 
nal fields. The latter perspective is taken in Ref. [2y| 
for the temperature field. We show in this section that 
the former is equivalent and more closely resembles the 
derivation required for polarization in 21 

For notational simplicity, we will derive the contribu- 
tion to the spectrum due only to the modulated piece of 
the model of Eq. ([T]) so that 



$(x) = g(x) [1 + /i(x)] . 



(Al) 



The unmodulated piece can be added to these results 
without loss of generality. 

The Fourier moments of $ are given by a convolution 
of the two fields g and h, 

m ^ g{k) + J ^g(k')Mk - k') , (A2) 



and the two point function of the field $ becomes 
(a>*(k)$(k')) ^i27rf6(k-k')Pgik) 

+ [Pg{k)+Pgik')]h{k' -k) 

(Pk 



(A3) 



-Pg{k)h*{k-k)h{k! -k) 



Under the Sachs- Wolfe approximation of Eq. ([5]), the 
spherical harmonic decomposition of the temperature 
field yields 



(Pk 

(2^ 



$(k)47rz^j,(fci?„c)r;„(k), (A4) 



which implies a two-point function of the temperature 
harmonics of 



I'm'} — t: 



(fik f (fk 



(27r)3 J (27r) 



($*(k)$(k')> 



x{ATifi^'-'ji{kD,,,)j,,{k'D,,,) 

xr,:„(k)rr™Kk')- 



(A5) 



Eqs. (jA3[) and (|A5p provide the analytic framework for 
modulated temperature two-point statistics. 



To make these considerations concrete, let us explicitly 
evaluate Eq. (|A5p for dipole modulation. Dipole mod- 
ulation can be taken to be a local approximation to a 
superhorizon scale Fourier mode. 



3 1 

/i(x) = wiW— — — — sin(ko-x) 



An fcoDroc 

fY z 

47r Dy^r. ' 



(A6) 



where we choose ko || z. The Fourier representation of 
the field becomes 

Now let us substitute this expression into Eq. (|A3p . 
The first term is zcroth order in h and is just the standard 
angular power spectrum 



(A8) 



where 



C, 



9 1 2^ 



(4^)j2(fcA-oc). 



The second term is first order in h and is given by 
1 f (Pk 



- 



9 J (27r)3 



d^fc' [Pg{k) + Pg{k')] 



2i V 47r fcoi^roc 
x[,5(k'-k-ko) -(5(k'-k + ko)] 

xi4nfi''-'jeikD,,,)je>ik'Drcc). 



The simplification of these expressions requires 
j(/{k'D^cc)yi'm'{^') to be expressed in terms of k where 
k' = k ± ko. Assuming that ko/k <^ 1, we can expand, 



1 



koD V 47r 
1 /' 3 



3 . /I, ,1 T^^^r / k + ako , 



|k-hako| 



koD V 47r 



n{kD)Yi^,{k) 



~ (^Rem'^^je+i{kD)Yi+i^,n{k) 
+ a^lm"'jf-i(fc^)^^-i,™(k), 
where the coupling matrix R\1^^ is defined as [3£ 

-{-ly 



(A9) 



(2£ + l)(2^i + l)(24 + l) 
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With these relations, the first order term in h becomes 



/rp* \(1) r nl.f ^TT I nl.f'^TT 



(All) 
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By virtue of the Wigner 3j symbol in the coupling matrix, 
the dipole modulation couples i to £±1 with a strength 
that is linear in h or wi. 

Likewise, the same operations on the term quadratic 
in h yield 



(A12) 



The quadratic term in h or wi couples £ to £±2. These 
results are the same as obtained in Rcf. 12611. 



2. Modulated Polarization Field 

For predicting polarization statistics, the direct angu- 
lar modulation approach of p6j does not apply given that 
the sources of the polarization field are modulated and 
not the field itself. 

The quadrupole source comes from the temperature 
field through Eq. Q 



(A13) 



where recall x' = x + ADh' (see Fig. [5]). 
Decomposing $ into Fourier modes yields 

^-^4^j2(fcAI?)r;,„(k)e^'^-. (A14) 

The Stokes parameters of the polarization are given by 
Eq. dED 



iQ±iU)ih) (A15) 
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m 

10 J dD 3 J (27r)3 ^ ' 

£' 

X ^ y2™.(k)r;„,(k)r,,™,(n)±2y2™.(n). 

7n2m' 

Recoupling the last line of Eq. (|A15p . we have 
(2f + 1)5 fe' 2 £ \ ( t' 2 t 



47r V t2 ±2 M 

x±2r,™(n)y,:„(k) 



Decomposing into the spin-2 spherical harmonics, 
Ei„. ± iBim. = - y dn±2r;„(n)(Q ± zf/)(n) 



10 y " Sy (27r)3 

Y.{\T,fn(kAD)Y,%''],,(kD) (A17) 



47r V =F2 ±2 / I 



Note that the i?-modcs vanish since the second Wigner 
3j symbol is zero if (£ + £' + 2) is odd and 



2 -2 j I -2 2 



(A18) 



This is to be expected since all scalar sources, even prod- 
ucts of fundamental scalars, generate only i5-modes. 
Paralleling the modulated temperature derivation of 
11 we find that there is the zeroth order piece coming 
from the unmodulated field that reproduces the standard 
result i38|]. 



EE 



where 



with 



I,(t) = J dD^e-'MkAD)ef'(kD), 

and the polarization radial function 



l3i£ + 2y.Mx) 



(A19) 



(A20) 



(A21) 



(A22) 



y 8(^-2)! x2 " 
The linear piece in h of the dipole modulation becomes 



/p* p \(i) - (47^)"^! 



1 



6 2 V 47r koD^oc 



(27r)3 J (27r)3 
X (27r)3 [5(k' - k - ko) - (5(k' - k + ko)] 



dr 



dD — e--j2{kAD)ef\kD) 

'-^-ir,™(k)yA„,(k')- 



(A23) 



(A16) 



Assuming, as in the temperature case, that fc/feg <C 1, 
we can again Taylor expand the spherical Bessel functions 
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and the spherical harmonic which, after some algebra, 
reduces Eq. (|A23p to 



where 



r,l.« 



I'm 



'(f + 2)!(^-2)! 
(^'-2)!(^ + 2)!' 



{5t,i 



'+1 



i)Bt 



A, 

■n, 

(A24) 



B, 



and 



Mk) 
Ke{k) 



dD—e 
dD 



dk k^Pg 



. D 
AD 



ATTh{k)J,{k) , 



ATih{k)Kt{k), (A25) 



■n{kAD)ef\kD) 



dD^e-^^n{kAD)ef'^ {kD) , (A26) 



with Ii{k) defined in Eq. (|A21|) . Like the temperature 
field, the coupling of ^ to £±1 is linear in h. The quadratic 
term can be similarly calculated but is more cumbersome 
and we omit the derivation here. 

For reasonable reionization histories, the dominant 
term in Eq. (|A24[) is the first term. Comparing Eq. (|A20p 
with ((A26)) we see that to order of magnitude. 



A, 



D 



Dr 



-c, 



EE 



(A27) 



where Z? is a typical distance to reionization. Com- 
pared with the modulation of the temperature spectrum 
in Eq. (jAll[) . the modulation of the polarization spec- 
trum is suppressed by a factor of D / D^^c- This factor is 
due to the lower level of spatial modulation at typical dis- 
tances to reionization in the dipolar model of Eq. (|A6p . 
Typically, it has a value of D/D^ac ~ 2/3. 

For the TE two point function, the same operations 
yield 



{Thrift 




(€'-2)!(^ + 2)! 
i'-i)Lui] -I- Cj,^ 



(A28) 



with 



Hf, = 





r dk k^Pg 


%/54 j 


k 27r2 




r dk k^Pg 


V54 j 


k 27r2 



ji{kDrcc)^TTJeik) , 



jeikD,,,)ATrKe,{k), (A29) 



where Cf^ is the usual unmodulated cross power spec- 
trum 



TE 



dk k^Pgik) 



/54 J k 27r2 
and It{k) is defined in Eq. (|A2ip 



47rj,(fcAec)/Kfc) , (A30) 



3. Comparison with Simulations 

In this section, we test the 3D simulation method 
against the analytic results of !jA II and lA 21 The simula- 
tions are performed according to the algorithm described 
in Section III Dl with a box of size that is four times the 
diameter of the last scattering surface and a spatial res- 
olution defined by the 512^ pixels. 

We first test the simulations against the usual results 
given by Eq. ((X9|) . (|A20|) . and (|A30|) for an unmodu- 
lated scale invariant potential spectrum. The temper- 
ature power spectra averaged over multiple realizations 
agrees with the analytic results within the sample vari- 
ance 



TT 



(jTT 



{2^ + l)N 



(A31) 



where N is the number of simulated skies. 

In order to predict the i?-modes, we calculate the 
quadrupolar temperature sources for 80 radii along each 
line of sight, and in every unit vector n corresponding to 
12, 288 HEALPIX pixel centers at resolution A^^^ido = 32. 
The small evolution of perturbations of the potential field 
between recombination and reionization is accounted for 
via the transfer function. 

We compute the E modes for different reionization his- 
tories with ionization fraction equal to 1 out to a redshift 
where the optical depth r = {0.0279, 0.055, 0.089} and 
zero beyond this redshift. In all cases, the distribution 
of numerical spectra are consistent with cosmic variance 
for i < 10. For larger multipoles there is a diminution 
of power that is consistent with the Nyquist frequency of 
the box, but this does not affect any of our results. 

Finally, we check the dipole modulation model against 
the analytic results of the previous sections. We showed 
there that a dipole modulation gives rise to couplings be- 
tween £ and £' = ^ ± 1 modes in the TT, EE, TE and 
ET two point functions. As an example of these corre- 
lations, we compare simulation results with the analytic 
calculation in Fig. 1151 We again find agreement within 
sample variance expectations. 



4. Dipolar Modulation Estimators 



From the analytic calculations of i|A II and lA 21 we can 
construct unbiased estimators of the dipole modulation 
amplitude wi. In the presence of a dipole modulation, 
the expectation values of temperature and polarization 
modes £ and £' = ^ ± 1 are correlated. Each pair of 
multipoles then forms an estimator of wi 



fXY r/U 



(A32) 



where XY G {TT,EE,ET,TE}, Rf+i^^ are coupling 
matrices defined in Eq. (jAlOp and the coefficients 
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FIG. 15: Comparison between analytical solution (dotted 
line) and 545 3D simulations (solid line) of the dipole modula- 
tion model with wi = 0.2. The two point correlation between 
the quadrupole and octopole m moments are shown for the 
X and Y fields: XY = TT, EE, TE, ET. Error bars reflect 
the sample variance expected for 545 skies. 



are given by 



(A33) 



such that (w^J^) = wi- Specifically, Eqs. ([XTT]) . (|X24)) . 
and (|X28|) give 



j-TT _r<TT 



riTT 



fEE 



cTE 



ET 
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(A34) 
1), 



TE 



Note that while individual wfj^^ estimators are complex, 
the joint estimator is real if -|-m and —m are combined 
with equal weights. 

We employ weights that minimize the variance of the 
joint estimator. Given a model with wi = 0, the variance 
is 



Here c'^xy x'v ^'"^ components of the matrix C given 

by 



XY,X'Y' Im 



XY aX'Y' '^XY,X'Y' 



(A37) 



■'XY,X'Y' 



fXY fX'Y' 



(A38) 



where Cf^ are the standard angular power spectra co- 
efficients. 

Imposing the normalization constraint on the weights, 
given by Eq. (jA36p . we minimize the variance with 



^XY^X'Y' X^fmP'^^lxyjf'y 



, (A39) 



which is the usual inverse covariance weight of covarying 
estimators. Likewise in the case that just one of the XY 
combinations is used this general form reduces to 



-XY ^"^ Cf^Cf/i 



( fXY nil 



(A40) 



which is the usual inverse variance weighted estimator. 
Note that the sum over (. can be restricted to I < ^max 
for modulation models which only affect large scales. 

Finally, given instrumental noise power spectra Ng^^ , 
the variance and minimum variance weights are modified 
by including the spectrum into the covariance matrix 



where Ag, Bg_i+i, Hg and Li^i+i are given by Eqs. (|A25[) 
and (|X2g|) . 

Each estimator carries a large sample variance and so 
one can combine them to form the joint estimator 



EE^ 



XY -XY 



(A35) 



XY 



where Af^ are the weights assigned to each one. To 
maintain the condition that the estimator is unbiased 
{wij ~ u'l, the weights must satisfy 



■'XY,X'Y' — fXY fX'Y' 

J I J I 



(A41) 



where we typically assume that the noise cross power 
spectrum TE vanishes. 



E E ^^'i' ~ ^ ■ 



(A36) 



XY hn 
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